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Abstract. We present recent QCD results from the D0 experiment at the Fermilab Tevatron Col- 
lider in pp collisions at = 1.96TeV. Results are presented for the inclusive jet and dijet cross 
sections, a measurement of dijet azimuthal decorrelations, studies of elastic scattering, and a search 
for diffractively produced Z bosons. 
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l/"") . The center-of-mass energy for proton-antiproton {pp) collisions was increased from 

y/s = 1.8TeV to 1.96 TeV in Run II of the Tevatron Collider at Fermi National Accel- 
erator Laboratory. This article presents QCD results measured with the D0 detector [1], 
based on data sets corresponding to integrated luminosities of ~ 150pb _1 . These are 
measurements of the inclusive jet and the dijet cross section, studies of azimuthal decor- 
relations in dijet production, a search for diffractively produced Z bosons, and a mea- 
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surement of the dN/d\t\ distribution in elastic pp scattering. 



The production rates of particle jets with large transverse momentum (px) provide 
sensitive tests for the predictions of perturbative QCD (pQCD) and give information 
on the non-perturbative structure of the proton as parameterized in the parton density 
functions (PDFs). The increased center-of-mass energy of y/s = 1.96 TeV in Run II 
results in an increased jet cross section, as compared to y/s = 1.8 TeV in Run I. Fig. 1 
shows the ratio of the pQCD predictions for the inclusive jet cross section for both 
center-of-mass energies. At pj above 400 GeV the cross section increases by more than 
a factor of two. At high integrated luminosity this increases the accessible pj range 
and extends the discovery potential for new physics phenomena. The increased cross 
section at high pj will also allow to test the dynamics of the strong interaction in 
unprecedented regions and constrain the PDFs at highest momentum fractions x. The 
sensitivity of the jet cross section to the PDFs can be seen from the decomposition of the 
total jet cross section into partonic subprocesses. This is shown in Fig. 1 (right) for the 
inclusive jet cross section as a function of pr- Up to pj = 200 GeV (where the PDFs are 
probed at x ~ 0.4), the largest fraction of the jet cross section comes from gluon-induced 
subprocesses. There are still significant contributions, even at high pj, from quark-gluon 
scattering, giving direct sensitivity to the gluon density in the proton even at highest x. 




FIGURE 1. The ratio for the NLO pQCD predictions for the central inclusive jet cross section at 
different center-of-mass energies as a function of jet pj (left). The decomposition of the central inclusive 
jet cross section into different partonic subprocesses as a function of jet pj (right). 



A major difference between the Tevatron jet measurements in Run I and Run II is 
the jet definition: the cone algorithm used in Run I was not infrared safe and reliable 
pQCD predictions were not possible [2]. The "Run II cone algorithm" [3] is still an 
iterative seed-based cone algorithm, but it uses midpoints between pairs of proto-jets 
as additional seeds which makes the procedure infrared safe [2]. A further difference 
between both jet algorithms is the chosen "recombination scheme". The Run I algorithm 
used the "£y-scheme" (or "Snowmass convention" [4]) in which the transverse jet 
energy Ej was defined as the scalar sum of the particle Ej's and the invariant jet mass 
was defined to be zero. Consequently, in Run I, the jet kinematics were described by 
Ej and pseudorapidity, r\ = — ln(tan(0/2)). In Run II particles are combined into jets 
using the "E-scheme" (addition of four-vectors). The it-scheme results in massive jets 
which are described by pj and rapidity, y = \ In ( (E + p z ) / (E — p z ) ) . The E-scheme has 
the advantage that the kinematic boundary pf 3 * = y/s/2 is stable, independent of the 
number of particles in a jet (which was not the case in the £V-scheme). This is desirable 
for resummed calculations [5]. 

All jet measurements presented here are using a data sample corresponding to an 
integrated luminosity of ~ 150pb The cone radius (in y and space) was chosen to 
be i?cone = 0.7. Jet detection was based on a compensating, finely segmented, liquid- 
argon and uranium calorimeter that provided nearly full solid-angle coverage. Events 
were acquired using multiple-stage inclusive-jet triggers. Analysis regions were defined 
based on the jet with largest pj in an event (p™ ax ) with the requirement that the 
trigger efficiency be at least 99%. The position of the pp interaction was reconstructed 
using a tracking system consisting of silicon microstrip detectors and scintillating fibers 
located within a 2T solenoidal magnet. The pj of each jet was corrected for calorimeter 
showering effects, overlaps due to multiple interactions and event pile-up, calorimeter 
noise effects, and the energy response of the calorimeter. The data were corrected for 
selection efficiencies and for migrations due to resolution. The largest experimental 
uncertainty is due to the jet energy calibration. 
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FIGURE 2. The inclusive jet cross section as a function of pr in different rapidity regions (a). Only 
statistical errors are shown. Figures (b), (c), (d) show the ratios of data and the NLO pQCD predictions. 
The error bars indicate the statistical errors and the total experimental uncertainty is displayed by the lines. 
Theoretical uncertainties due to the PDFs are shown as the shaded bands. 



Inclusive Jet and Dijet Jet Cross Sections 

The inclusive jet cross section is measured as a function of pj and y. Fig. 2 (a) 
shows the pj dependence of the inclusive jet cross section in three rapidity regions 
in the range < \y\ < 2.4 (with statistical errors only). The cross section is falling 
by more than six orders of magnitude between pj = lOOGeV and pj = 500 GeV. At 
forward rapidities the pr dependence is much stronger. Next-to-leading order (NLO) 
pQCD predictions are overlaid on the data. The NLO calculations [6] were computed 
for renormalization and factorization scales \x r = \lf = 0.5 /?™ ax using the CTEQ6M [7] 
PDFs and a s (Mz) = 0.1 18. The maximum distance of particles within a jet was limited 
to AR < R sep ■ i?cone with i? sep = 1.3 [8]. The ratio of data over theory is shown in Figs. 2 
(b), (c), (d). Uncertainties in the NLO calculations due to the PDFs are indicated by the 
grey bands; the experimental uncertainties are displayed as lines. The latter increase with 
Pt, especially at large rapidities. Theory has good agreement with data given the large 
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FIGURE 3. The dijet cross section as a function of My at central rapidities. In Fig. (a) the data are 
shown with the total experimental uncertainties. The statistical uncertainties are displayed by the inner 
error bars. Fig. (b) shows the ratio of data and the NLO pQCD predictions. The error bars indicate the 
statistical errors and the total experimental uncertainty is displayed by the lines. Theoretical uncertainties 
due to the PDFs are shown as the shaded bands. 



uncertainties. The dijet cross section at central rapidities (\y\ < 0.5) is shown in Fig. 3 (a) 
as a function of the invariant dijet mass My in the range 150 < My < 1400 GeV. The ratio 
of data and theory is displayed in Fig. 3 (b). Within the large experimental uncertainties 
theory gives a good description of the data. 

It is clear that the current jet cross section measurements with their large uncertainties 
can not constrain the PDFs. Improvements in the jet energy calibration and a reduction 
of the corresponding uncertainty are still needed. It has to be noted that at the time this 
article is written significant improvements have already been achieved and further work 
is in progress. 



Dijet Azimuthal Decor relations 

As long as the uncertainties of the jet energy calibration are still large it is convenient 
to investigate observables which have a smaller dependence on the jet energy calibration 
but which are still sensitive to QCD effects. The decorrelation of the azimuthal angles 
of the two leading pj jets, A0 dijet = l^jeti — 0jet2|> lS suc h an observable. In lowest 
order pQCD both jets in a dijet event have equal pj and correlated azimuthal angles 
dijet = Additional soft radiation causes small azimuthal decorrelations, whereas 
dijet significantly lower than % is evidence of additional hard radiation with high 
Pt- Exclusive three-jet production populates 2tz/3 < A0(jij et < % while smaller values 
of A0 dijet require additional radiation such as a fourth jet in an event. Distributions 
in A0 dijet provide an ideal testing ground for higher-order pQCD predictions without 
requiring the reconstruction of additional jets and offer a way to examine the transition 
between soft and hard QCD processes based on a single observable. 

The observable is defined as the differential dijet cross section in A0dijet> nor- 
malized by the inclusive dijet cross section integrated over the same phase space 
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FIGURE 4. Left: The A<j> 

dijet distributions in four regions of /?™ ax . The solid (dashed) lines show the 
NLO (LO) pQCD predictions. Right: Ratios of data to the NLO pQCD calculation for different regions 
of pf ax . Theoretical uncertainties due to variation of \i r and \ij are shown as the shaded regions; the 
uncertainty due to the PDFs is indicated by the solid lines. The points at large A(j> <jij e t are excluded because 
the calculation has non-physical behavior near the divergence at n. 



(l/°dijet) (<iOdijet/^A0 dijet)- The corrected data [9] are presented in Fig. 4 (left) as a 
function of A0 dijet in f° ur ranges of p™ ax . The inner error bars represent the statistical 
uncertainties and the outer error bars correspond to the total experimental uncertainties. 
The spectra are strongly peaked at A0 dijet ~ ft and the peaks are narrower at larger values 
of 

pmax rp^ e ^ ata - n ^ (left) are compared to the pQCD predictions in leading order 
(LO) and NLO [10] for \i r = jlf = 0.5p^ ax using CTEQ6.1M [7] PDFs. The (N)LO 
predictions are computed as the ratio of the (N)LO predictions for 2 — ► 3 processes 
and 2^2 processes. Due to kinematic constraints, the LO calculation is restricted to 
A0 dijet > 2^/3. At A0 dijet = ft where the third jet becomes soft the LO predictions di- 
verge. NLO pQCD provides a good description of the data. This is emphasized in Fig. 4 
(right) which displays the ratio of data and NLO. There is good agreement within a 
small offset of 5-10% relative to unity. Theoretical uncertainties due to the PDFs [7] are 
estimated to be below 20%. Also shown are the effects of \x r and \Xf variations in the 
range 0.25 p™ ax < \x r j < pf AX . The large scale dependence for A0 dijet < 2^/3 occurs 
because the NLO calculation only receives contributions from tree-level four-parton fi- 
nal states in this regime. NLO pQCD fails to describe the data in the region A0 dijet ~ ft 
which is dominated by soft processes and therefore can not be described in fixed- order 
perturbation theory. pQCD results at large A^) dijet m Figs. 4 were excluded because the 
calculations have non-physical behavior near the divergence at n and were not stable. 
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FIGURE 5. The A<p di j et distributions in different /?™ ax ranges. Results from HERWIG and PYTHIA are 
overlaid on the data. The left Figure displays the full A<p dijet range and the right Figure emphasizes the 
region close to %. 



Monte Carlo event generators, such as HERWIG and pythia, use 2^2 LO pQCD 
matrix elements with phenomenological parton-shower models to simulate higher order 
QCD effects. HERWIG and PYTHIA results for default parameters, CTEQ6L [7] PDFs, 
and a s (Mz) =0.118, are compared to the data in Fig. 5 over the whole A0 dijet range 
(left) and in greater detail in the region A(j) di j et > 13 ft/ 16 (right). The default version of 
HERWIG gives a good description of the data over the whole A0 dijet range in all p™ ax 
regions. It is slightly below the data around A(f) di j et ps 7n/8 and slightly narrower peaked 
at Tt. The default version of PYTHIA does not describe the data. The distribution is too 
narrowly peaked at A(j) dijet ~ n and lies significantly below the data over most of the 
A0 dijet range in all pf ax regions. 

To investigate the possibilities of tuning PYTHIA we focus on the impact of the initial- 
state radiation (ISR) parton shower and the parameters by which the ISR shower can 
be adjusted. The maximum allowed pj in the ISR shower is controlled by the product 
of the hard scattering scale (p^) and the parameter PARP ( 67 ) . Increasing this cut-off 
from the current default of 1 .0 to the earlier default of 4.0 leads to significant changes of 
the PYTHIA prediction for A(j) dijet (the shaded bands in Fig. 5, left, and the dotted line in 
Fig. 5, right). At low A0 dijet the description of the data becomes very good but at large 
A0 dijet this parameter has not enough effect to bring PYTHIA close to the data. The best 
description of the low A0 dijet region is obtained for PARP ( 67 ) =2 . 5 (not shown). 

Further ISR-related parameters are tested if they have impact on the A0 dijet distri- 
bution. These are the scale factor, Xn, for the renormalization scale for <x s in the ISR 




FIGURE 6. Left: The A0 di j et distributions in different /?™ ax ranges. Results from PYTHIA are overlaid 
on the data for different parameter settings of the ISR shower. Right: Non-perturbative corrections for the 
A<t> 

dijet distribution. Shown are the hadronization corrections (solid line) and the effects from underlying 
event as determined using PYTHIA. 



shower, PARP ( 64 ) , and the primordial kj of partons in the proton: the central value 
of the gaussian distribution, PARP (91), and the upper cut-off, PARP (93). We have 
lowered the factor for the renormalization scale to PARP ( 64 ) = . 5 (D=1.0) which 
increases the value of a s according to the RGE. We have alternatively increased the pri- 
mordial k T from 1 GeV to 4 GeV, PARP ( 9 1 ) =4 . (D=l . ) , and the upper cut-off 
of the gaussian distribution from 5 GeV to 8 GeV, PARP ( 93 ) =8 . (D=5 . ) . These 
parameter variations have no effect on the region at low dijet and the effects at large 
A(j> Hijet are shown in Fig. 6 (left). Both effects are very small: while the scale factor has 
almost zero influence (dashed-dotted line) there is some small change for the increased 
primordial kj (solid line) which manifests itself, however, only very close to the peak 
region and only at lower p™ ax . 

Non-perturbative corrections for the A0 dijet distribution are investigated in Fig. 6 
(right) using PYTHIA. Hadronization corrections are defined as the ratio of the observ- 
able, defined on the level of stable particles and on the level of partons. The underlying 
event correction is defined as the ratio of the predictions including underlying event and 
with the underlying event switched off. Both corrections are below 5%, decreasing with 
pmax rp^g sma i mess f these corrections allows to directly compare the measurement to 
purely perturbative QCD predictions and it also allows to interpret the discrepancies of 
the Monte Carlo event generators in terms of perturbative processes. 
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FIGURE 7. The logarithm of the reconstructed energy in the direction of the outgoing antiproton in the 
range — 5.3 < tj < —2.6. Compared are the distributions for events with no visible interactions (line) and 
events in which the proton and the anti-proton dissociate (full and open markers). The areas are normalised 
to unity. Rapidity gap candidates are selected by applying a cut at 10 GeV. 



So far we have discussed inelastic processes with large transverse momentum. However, 
only a small contribution to the total pp cross section comes from such hard inelastic 
processes. About 40% of the total pp cross section is elastic scattering and diffraction. 
In elastic scattering both protons are scattered with no momentum loss under a small 
angle and emerge intact. No further particles are produced. In single diffraction only 
one of the protons remains intact and the other proton dissociates. While additional 
particles can be produced, there is often a large rapidity interval towards the intact proton 
in which no particles are emitted (rapidity gap). Both elastic scattering and diffractive 
processes can be identified experimentally by either measuring the scattered protons or 
by reconstructing the rapidity gap as an angular region with very little energy in the 
calorimeter. In the following sections we present a search for diffractively produced 
Z bosons in the muon decay channel and a preliminary spectrum of events in elastic 
scattering as a function of the four-momentum transfer to the proton. 



D0 collected an event sample of about 21 pb 1 during Run I suitable for diffractive 
physics studies. Nine diffractively produced Z — > e + e~ candidates were observed in this 
data sample [13]. The current search for Z bosons, produced in single diffraction, uses 
a dataset of ~ 110 pb . The analysis is based on the experimental reconstruction of a 
rapidity gap close to the direction of the final-state (anti-) proton and the reconstruction 
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FIGURE 8. The invariant di-muon invariant mass distribution for Z boson candidates with a single 
rapidity gap (left) and with no rapidity gap (right). For details on the rapidity gap definition, please see 
the text. 



of the Z via its decay into two muons with opposite charge. Both muons were required 
to have pj > 15 GeV. One of the muons must be isolated in the calorimeter and in the 
central tracking detector. Muons from cosmic rays were suppressed by the requirement 
that their tracks stem from the vertex of the pp interaction and that their acolinearity 
(cj> = A0 MjU + A0 M|U - In) is larger than 0.05 [14, 15]. The search of the rapidity gaps 
uses two scintillating detectors, one on each side of the interaction region, which cover 
the pseudorapidity range 2.7 < 1 77 1 < 4.4. Calorimeter energies towards the directions of 
the outgoing protons (in the range 2.6 < |rj| < 5.3) are summed separately on each side 
with cell thresholds of 100 MeV in the electromagnetic layers and 200 MeV in the fine 
hadronic layer. The distribution of the logarithm of the energy sum is shown in Fig. 7 
for bunch crossings with no visible interactions. Based on a random trigger these events 
are expected to have similar topologies as rapidity gap events with no activity in the 
direction of the outgoing antiproton. Also shown is the energy distribution for a sample 
of minimum bias events. A third event sample is defined by the requirement of (at least) 
one jet with pj > 25 GeV. Events are excluded when the leading pj jet is in the region 
1 7] I > 2.4. Both the minimum bias and the jet samples are dominated by non-diffractive 
events in which both protons dissociate. Events where the antiproton remains intact and 
events with antiproton dissociation can be separated by requiring the energy sum to be 
below 10 GeV. The final selection of single diffractive Z boson candidates requires that 
on one side the scintillator is off and the energy sum is below 10 GeV while on the other 
side the scintillator is on and the energy sum is larger than 10 GeV. 

The invariant di-muon mass distribution is shown in Fig. 8 for Z boson candidates 
in events with a rapidity gap (left) and in events that fail the rapidity gap selection 
(right). The latter are strong candidates for non-diffractively produced Z bosons. Their 
mass distribution has a strong peak at the Z mass and small background. The di-muon 
mass distribution for the candidates for diffractively produced Z bosons has very similar 
shape. Additional studies on the purity and the background for the rapidity gap selection 
are, however, still needed for quantitative statements and to establish the evidence for a 
diffractive Z — > signal. 
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FIGURE 9. The differential elastic pp scattering cross section as a function of the four- momentum 
transfer squared \t\. Measurements from experiments at different center-of-mass energies are shown 
(datasets at different energies are scaled for visual separation). The preliminary D0 data are arbitrarily 
normalized for comparison of the slope with the E710 data at y/s = 1 .8 TeV. Theoretical predictions from 
Block et al. are overlaid on the data. 



Elastic Scattering 

In the D0 experiment an intact outgoing (anti-) proton can be detected in the Forward 
Proton Detector (FPD) [16], which consists of position detectors along the beam line 
in conjunction with accelerator magnets. In total, the FPD has nine spectrometers, each 
comprising two scintillating fiber tracking detectors that can be moved to within a few 
millimeters of the beam. The kinematics of the outgoing (anti-) proton is described by 
the four-momentum transfer squared, t and the fractional longitudinal momentum loss 
^ . Both variables \t\ and t, can be computed from the track four- vectors reconstructed in 
the FPD. 

We have measured the \t\ distribution in elastic scattering using the FPD. Elastic 
events were triggered by two spectrometers which are above (below) the beam on the 
side of the outgoing antiprotons (protons). The spectrometer on the proton side was used 
to measure the track position. The trigger uses time-of-flight information, and vetoes 
events with hits in the scintillators covering 5.2 < |rj| < 5.9. The spectrometers can 
reconstruct protons in the kinematic range 0.8 < \t\ < 4.0 GeV 2 , depending on detector 
position and accelerator conditions. The preliminary results are presented in the range 
0.92 < \t\ < 1.34 GeV 2 . 

A Monte Carlo event generator was used to propagate the outgoing (anti-) protons in 



elastic scattering through the magnetic field of the Tevatron and to simulate the response 
of the FPD based on the expected detector resolution [17]. The geometric acceptance 
of the FPD was modeled by a similar event generator [18] as a function \t\. Remaining 
backgrounds from beam halo [19] and from single diffraction are subtracted separately. 

The preliminary dN/d\t\ distribution for elastic pp scattering is shown in Fig. 9. Mea- 
surements of the differential cross section do/d\t\ at different center-of-mass energies 
from the ISR, UA4, and E710 experiments are also shown [20]. The normalization of 
the D0 data is arbitrary, but it is visible that the shape of the D0 data in the range 
0.8 < \t\ < 1.4GeV is very different from the shape of the E710 data measured at 
y/s = 1.8 TeV and at \t\ < 0.65 GeV. A theoretical prediction by Block et al. [21] is 
compared to the data at highest center-of-mass energies. The model predicts a transition 
between the different \t\ regions and give a good description of the data over the whole 
\t\ range. 
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